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ABSTRACT 

The Planck satellite provides a set of all-sky maps spanning nine frequencies from 30 GHz to 857 GHz, with an unprecedented combination of 
sensitivity and angular resolution. Planets, minor bodies, and diffuse interplanetary dust (IPD) contribute to the submillimetre and millimetre sky 
emission it observed. The diffuse emission can be effectively separated from Galactic and other emissions, because Planck views a given point 
on the distant celestial sphere multiple times, through different columns of IPD. We use the Planck data to investigate the behaviour of Zodiacal 
emission over the whole sky in the sub-millimetre and millimetre. We fit the COBE Zodiacal model to the Planck data to find the emissivities of the 
various components of this model - a diffuse cloud, three asteroidal dust bands, a circumsolar ring, and an Earth-trailing feature. The emissivity of 
the diffuse Zodiacal cloud decreases with increasing wavelength, as expected from earlier analyses. The emissivities of the dust bands, however, 
decrease less rapidly, indicating that the properties of the grains in the bands are different than those in the diffuse cloud. As part of the analysis, 
we fit the small amount of Galactic emission seen through the instrument's far sidelobes and place limits on possible contamination of the CMB 
results from both Zodiacal emission and Galactic emission seen through these far sidelobes. When necessary, these results are used in the Planck 
pipeline to make maps with Zodiacal emission and far sidelobes removed. We show that the spectrum of the Zodiacal correction to the CMB maps 
is small compared to the Planck CMB temperature power spectrum. 

Key words. Zodiacal Dust - Interplanetary Medium ~ Infrared: Solar System - Submillimetre - (Cosmology): Cosmic Microwave Background 



1. Introduction 

This paper, one of a set associated with the 2013 release of 
data from the Planck^ mission (Planck Collaboration I 2011), 



* Corresponding author: Ken Ganga - Ken . GangaOapc . 
univ-paris-diderot . fr 

' Planck (http://www.esa.int/Planck) is a project of the European 
Space Agency (ESA) with instruments provided by two scientific con- 



describes the measurement of Zodiacal emission using Planck 
data. 

Zodiacal light, or the reflection of sunlight from small dust 
particles in our Solar System, can be seen by eye at dawn or dusk 



sortia funded by ESA member states (in particular the lead countries 
France and Italy), with contributions from NASA (USA) and telescope 
reflectors provided by a collaboration between ESA and a scientific con- 
sortium led and funded by Denmark. 
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in dark locations, and contributes significantly to the diffuse sky 
brightness at optical and near-infrared wavelengths. The study 
of Zodiacal emission, or the thermal re-emission of absorbed 
energy from these interplanetary dust (IPD) particles has been 
enabled by the advent of infrared astronomical techniques, and 
it is now known to dominate the diffuse sky brightness over most 
of the sky between 10 and 50 yum (see, for example, Leinert et al. 
1997). 

Full-sky, infrared satellite surveys, in particular, have al- 
lowed us to begin to determine the structure of the density of 
the IPD (Hauser et al. 1984; Kelsall et al. 1998; Fixsen & Dwek 
2002; Pyo et al. 2010). One of the fall-sky models of Zodiacal 
emission most easily adapted for Planck and most commonly 
used at longer wavelengths is based on work from the Cosmic 
Background Explorer Diffuse Infrared Brightness Experiment 
(COB^/DIRBE) team (Kelsall et al. 1998, hereafter K98). 
Others are presented in Good et al. (1986), Rowan-Robinson 
et al. (1990, 1991), Jones & Rowan-Robinson (1993), Vrtilek 
& Hauser (1995), Wright (1998), and Rowan-Robinson & May 
(2012). The K98 model comprises the well-known diffuse cloud, 
three sets of dust bands first discovered by IRAS (Low et al. 
1984), and a circumsolar ring and Earth-trailing feature, hinted 
at in IRAS and confirmed in DIRBE (called a 'blob' in K98. See 
Reach et al. 1995, and references therein). 

Fixsen & Dwek (2002) have used data from the Far Infrared 
Absolute Spectrophotometer (FIRAS) to extend measurements 
of the diffuse cloud to longer wavelengths, but given its mod- 
est angular resolution and large uncertainties in the submillime- 
tre region, could not say more about the smaller angular-scale 
Zodiacal features. Planck's sensitivity allows it to detect and 
measure the emissivity of the diffuse Zodiacal cloud at long 
wavelengths, and its angular resolution also allows it to char- 
acterize the smaller-scale components of the Zodiacal. 

This paper continues as follows: in Sect. 2 we describe the 
relevant aspects of the Planck mission for this analysis, along 
with it's observation strategy and data processing; in Sect. 3 we 
describe how Planck detects Zodiacal emission; and in Sect. 4 
we briefly describe the COBE Zodiacal emission model. The fit 
to the Planck/HFl data is described in Sect. 5, and the results of 
the fit are discussed in Sect. 6. We conclude in Sect. 7. 



2. The Planck Mission 

Planck, comprised of the High Frequency Instrument, or HFI, 
and the Low Frequency Instrument, or LFI, was launched in 
May of 2009. The mission as a whole is described in Planck 
Collaboration I (2013). This work uses only data at frequencies 
of 100 GHz and higher. At these frequencies, Planck observed 
the entire sky in six broad frequency bands from 100 to 857 GHz, 
with corresponding angular resolutions from roughly 9.'7 to 
4.'6 (Planck CoUaboration VII 2013). 

2.1. Orbit, Scanning Strategy and Dates of Observation 

While Planck's orbit and scanning strategy are described in 
depth in Planck Collaboration I (201 1) and Planck Collaboration 
I (2013), we give a synopsis of the elements relevant to our anal- 
ysis here. 

Planck orbits around the second Sun-Earth Lagrange point, 
and is thus always close to the Ecliptic plane and about 1.01 AU 
from the Sun and 0.01 AU from the Earth. Its focal plane scans 
the sky once per minute, with each detector always observing 
on a circle approximately 85° from its spin axis. A simplified 



video showing this scanning strategy can be found at the ESA 
website^. In addition, and not visible in the video noted above, 
the spin axis traces the Sun-Earth vector, but with an additional 
"cycloid" component, so that in the Sun-Earth frame the spin 
axis is always 7.5° degrees from the Sun-Earth vector and cir- 
cles around it twice per year. This cycloid component results in 
differing total amounts of IPD in Planck's line of sight for differ- 
ent observations of the same point on the distant celestial sphere. 
This is shown schematically in Fig. I. 




Fig. 1. Schematic representation of the geometry of Planck's 
measurements, which shows that it can view different amounts 
of Zodiacal emission while looking at the same point on the dis- 
tant sky. The plane of the ecliptic is in the plane of the diagram. 
The Sun is in the center of the circles. The solid black line rep- 
resents the orbit of the Earth and Planck. The dashed line at the 
outer edge of the shaded ring represents the orbit of Jupiter, be- 
yond which we assume there is no contribution to the Zodiacal 
emission from IPD. Panel (a) shows a case where the phase of 
the scan cycloid and the location of the observed point on the 
sky yield two measurements for which the lines of sight through 
the IPD is roughly equal, and the same Zodiacal signal is seen. 
Panel (b) shows a case where the phase of the scan cycloid and 
the location of the observed point on the sky yield different to- 
tal columns of IPD along the lines of sight, and thus a different 
Zodiacal signal is seen in each of the two measurements. Note 
that this figure is highly stylized and not to scale. 

As nearly the entire sky is seen twice each year, the Planck 
team divides the observations into "surveys" of approximately 
six-month duration. The exact definition of the beginning and 
end of each survey was agreed upon within the Planck team. The 
basic characteristics are that each survey lasts about six months 
and covers a maximum of sky, with a minimum of overlap be- 
tween the beginnings and ends of the survey. During any single 
one of these surveys, some pixels near the ecliptic poles are ob- 
served multiple times, as are the pixels near the ecliptic plane 
which are seen both at the beginning and at the end of the sur- 
vey. The bulk of the sky, however, is observed only during well- 
defined periods, usually less than a single week. In Fig. 2 we 
show the Julian dates of observations of those pixels on the sky 
for which the observation times during survey 1 spanned one 
week or less. The analogous plot for survey 2 is similar in na- 



^ http : //sci . esa . int/ science- e/www/ob j ect/ index . cfm? 
fobjectid=47345. 
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ture, and the corresponding maps for surveys 3 and 4 are quite 
similar to those of surveys 1 and 2, as the scanning strategy for 
surveys 3 and 4 were almost identical to those of surveys 1 and 
2, respectively. 




2.45506e+06 Julian Date 2.455246+06 



Fig. 2. The Julian date of observation of pixels on the sky during 
survey 1, for a single detector, in Galactic coordinates. There are 
only very small differences between maps for different detectors. 
The grid lines show ecliptic coordinates, with the darker lines 
representing the ecliptic plane and the line of zero ecliptic longi- 
tude. Undefined pixels, which were either not observed at all, or 
which were observed multiple times over a period that spanned 
more than one week and are thus not used in this analysis, are 
shown as the uniform gray band. 



2.2. Data Processing 

The overall HFI data processing is described in Planck HFI 
Core Team (2011b) and Planck Collaboration VI (2013). Given 
the time-dependent nature of the Zodiacal signal seen with the 
Planck scanning strategy, this analysis is done using the individ- 
ual survey 1-4 maps. This allows us to exclude from the analysis 
regions of the sky and periods of time where the column of IPD 
viewed by Planck is not constant. 

The HFI instrument has a number of horns at each mea- 
surement frequency (Planck Collaboration VII 2013, Fig. 9). 
Working with individual horn maps, rather than the co-added fre- 
quency maps, allows us to adjust the response of each detector 
so that they are uniform for a source with a Zodiacal spectrum, 
rather than that used for a CMB spectrum, as is done in the stan- 
dard processing (Planck Collaboration IX 2013). At 100, 143, 
217, and 353 GHz, some of these horns contain two polariza- 
tion sensitive bolometers (PSBs; Jones et al. 2003). As we are 
not addressing polarisation here, for those horns that have PSBs, 
we combine the maps from each of the two detectors within the 
given horn with a simple average. 

As the evaluation of the model to be presented in Sect. 4 in- 
volves calculating emission from a number of points along each 
line of sight and summing them, the computations are time con- 
suming. To mitigate this to some extent, we use 13.'7 x 13.'7 
pixels, rather than the original 1 .'7 x 1 .'7 HFI pixels, reducing the 
number of map pixels from 50 million to a bit less than 800 thou- 
sand (that is, we use HEALPix pixels with A^side of 256 rather 
than 2048; Gorski et al. 2005). While this does reduce our sensi- 
tivity to finer scale structures, it is not a serious hindrance, since 
we will be making comparisons with DIRBE, which had a still 



larger beam. Smaller pixels will be used in future work as more 
detail is teased out of the data. 

Pre-launch estimates of Planck's ability to detect Zodiacal 
emission and an estimate of the possible level of contamination 
at the highest Planck frequencies were presented in Maris et al. 
(2006). More recent predictions have addressed the possibility of 
Zodiacal contamination at lower frequencies (Diego et al. 2010) 
and speculated that emission from dust in the outer Solar System 
might contribute to large-scale anomalies which have been re- 
ported in data from the Wilkinson Microwave Anisotropy Probe 
(WMAP) at large angular scales (Maris et al. 201 1; Hansen et al. 
2012). 



3. Detection 

The existence of the Zodiacal emission in the Planck maps 
is straightforward to demonstrate by exploiting the fact, noted 
above, that different Planck surveys often sample different 
columns of IPD while observing the same location on the distant 
celestial sphere. Fig. 3 shows the first (top) and second (middle) 
survey maps for the 857-1 detector, along with the difference 
of these two maps (bottom). Similar differences for all HFI fre- 
quencies, using data from all horns at each frequency, are shown 
in the left-hand panel of Fig. 10. There are at least three fea- 
tures that stand out in the map difference: (1) The scale has 
been reduced immensely, but the Galactic plane is still visible 
in this difference map; (2) the "arcs" at the top and the bottom 
of the difference map are the images of the Galactic centre as 
seen through the instrument's far sidelobes (FSLs); and (3) the 
Zodiacal emission can be seen as the variations following the 
ecliptic plane. We begin with the Zodiacal emission. 

Since these are difference maps, the Zodiacal emission is 
seen in the ecliptic plane as both positive and as negative, de- 
pending upon the relative geometry of the IPD and the Planck 
satellite when a given location was observed in the two surveys. 
While the amplitude of the Zodiacal emission signal is reduced 
in this differential process (this will be described in Sect. 4, and 
in particular Fig. 4; see also Maris et al. 2006), it has the ad- 
vantage that the bulk of the Galactic and extra-galactic signal, 
the main contaminants in the analysis at high frequencies, is ef- 
ficiently removed with little effort. What remains of the Galactic 
signal arises from effects such as beam asymmetries and imper- 
fections in the transfer function removal (Planck Collaboration 
VII 2013). To mitigate these residual Galactic effects, we gener- 
ally do not use data within 10 degrees of the Galactic plane in 
this analysis. We have confirmed that changing this cut to 5 or 
to 20 degrees does not change the conclusions of this analysis, 
and emphasize that contrary to other analyses, the fact that we 
are analysing difference maps makes this work less sensitive to 
Galactic contamination. 

It should be noted that the Zodiacal emission is much dim- 
mer than many other background components in the Planck data. 
Whereas the Zodiacal emission dominates the sky in some IRAS 
and COBE bands, this is never the case for Planck, with the cos- 
mic infrared background and Galactic dust dominating at high 
frequencies, and the CMB itself dominating at lower frequen- 
cies. This makes our differencing scheme appealing, but also re- 
stricts the analysis in some ways. It is, for example, difficult to 
look at individual scans, or slices of the sky, as has been done 
successfully with IRAS. We are obligated to use almost the en- 
tire sky, and use a model of the Zodiacal emission to interpret 
variations, rather than being able to directly interpret the total 
Zodiacal emission on the sky. 
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Fig. 3. Single-survey maps in Galactic coordinates for the 857-1 
detector. Top: Survey 1 map. Middle: Survey 2 map. Bottom: 
Survey 2 minus survey 1 difference map. This bottom image 
shows the Zodiacal emission and the residual Galactic emission 
effects discussed in the text. The units are MJy/sr, assuming a 
spectrum inversely proportional to frequency. Undefined pixels 
are shown in gray. These occur in pixels which either have not 
been observed during the survey, were observed during the pas- 
sage of a planet, or a small number of other events. In the top 
two plots, pixels which were observed over periods longer than a 
week were not masked, and thus the masked regions are smaller 
in the top two images than that in the difference (bottom) panel. 



4. Model 

The goal of this section is to describe the Zodiacal and far side- 
lobe templates we will create to fit to the map at the bottom of 
Fig. 3. 



4.1. Zodiacal Components 

The COBE/DIRBE Zodiacal emission model is described in 
depth in K98, but we review the salient parts here. 



4.1 .1 . Diffuse cloud 

The density of the diffuse IPD cloud, having both radial and ver- 
tical dependence, is taken to be of the form 



where 



Rc = ^jix-xof + (y-yof + iz- 

^ = \Zc\ /Rc, 



zof, 



= (x - xo) ■ sin (Qr) ■ sin (/«) 

- iy- yo) ■ cos (Qjj) ■ sin (/r) + {z- Zo) ■ cos(/r), 



(1) 

(2) 
(3) 

(4) 



and a, /3, y, /j, xq, yo, zo, ^r, and are parameters describing 
the location and shape of the cloud. This form (and others used 
elsewhere which are similar) is based on an approximation of 
a model of particles orbiting the Sun, accounting for drag from 
the Poynting-Robertson effect, and with a modified fan distribu- 
tion used to describe the changes in density above and below the 
plane of the ecliptic. See K98 for more details and references. 

The numerical values for the parameters can be found 
in K98, or from the LAMBDA website^^. This is shown, for both 
survey 1 and survey 2, as well as their difference, assuming an 
emissivity of one, in the bottom row of Fig. 4. 

4.1.2. Dust bands 

The Zodiacal dust bands were first seen by IRAS (Low et al. 
1984), and appear as pairs of bright, parallel bands equally 
spaced above and below the ecliptic plane. They were quickly 
associated with asteroid families, and then understood to be the 
relics of asteroid collisions or collapses (Dermott et al. 1984; 
Sykes & Greenberg 1986). Reach et al. (1997) study them in 
detail. 

The K98 model contains three bands'* called, reasonably 
enough, bands 1, 2 and 3. They appear at ecliptic latitudes 
around +1?4, +10°, and ±15°. IRAS, having higher angular res- 
olution than DIRBE, was able find more bands, dubbed a, y, 
E/F, G/H, J/K, and M/N (Sykes 1990), though some are more 
firnily detected than others. 

The K98 band 1, or IRAS band y, was originally associ- 
ated with the Eos family of asteroids (Dermott et al. 1984), 
but Grogan et al. (2001) called this into question and Nesvorny 
et al. (2003, hereafter N03) found better correspondence with the 
Veritas family of asteroids. 

K98 note that their band 2 is a blend of IRAS bands a 
and j6 (Sykes 1990). Sykes & Greenberg (1986) tentatively as- 
sociated the a band with the Themis family of asteroids, and 
Nesvorny et al. (2008, hereafter N08) has narrowed this associa- 
tion to a cluster within this family associated with the Beagle as- 
teroid. The /3 band was associated by Sykes & Greenberg (1986) 



^ Legacy Archive for Microwave Background Data Analysis - NASA 
(http : //lambda . gsf c. nasa . gov/). 

While three bands are described in the text of K98, there are actu- 
ally four in the code. There is, however, no ambiguity, as the density of 
the fourth band is set to zero in the code at the LAMBDA site. 
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Fig. 4. Estimated 857 GHz Survey Zodiacal emission templates based on the model of K98, in Galactic coordinates, for survey 1 
(left), survey 2 (center), and their difference (that is, survey 2 minus survey 1; on the right). From top to bottom we show dust band 
1, dust band 2, dust band 3, the circumsolar ring, the Earth-trailing feature, and finally the diffuse cloud. Note that the scales for the 
right-hand column are somewhat different from those of the left-hand and centre columns, and that the scale for the bottom row is 
different from that of the others. 
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with the Koronis family of asteroids, and N03 narrowed this to 
the Karin cluster within the Koronis family. Anticipating the dis- 
cussion in Sect. 6.4, we note also that the yS band appears brighter 
than the a band (Sykes 1988; Reach et al. 1997; Nesvorny et al. 
2008). 

K98 states that their band 3, that furthest from the ecliptic 
plane, has been associated with both the lo and Maria families 
of asteroids (Sykes 1988; Reach et al. 1997), corresponding to 
IRAS bands J/K and M/N. N03 have more recently noted, how- 
ever, that 4652 lannini and/or 845 Naema may be better asteroid 
associations for the J/K band-pair and that 1521 Seinajoki may 
work better for the M/N pair We summarize these associations 
in table 1 . 

For each of the three dust bands in the COBE model, the 
density is given by 

(R) . ^,-(^/*.)' ( 1 + (1 _ (5) 

where we have used a simplified notation based on that of K98, 
where one can also find the numerical values for the parameters. 
Note that Eqn. 5 matches the code used for the Zodiacal model 
(which can be found on the LAMBDA website), but that there 
is a factor of l/v^ difference between Eqn. 5 and Eqn. 8 of K98 
(vg is a shape parameter, not a frequency). Also, K98 assumed 
that the emissivities of the three sets of bands were all equal. We 
relax this assumption below and allow the emissivities of each of 
the sets of bands to be different. The bands are shown, assuming 
unit emissivity, as the first, second, and third rows in Fig. 4. 

4.1.3. Circumsolar ring and Eartli-trailing feature 

The functional form for the density of the circumsolar ring is 
taken to be 

riB (R) = nsR ■ e-<«-«"'>'/^'s«-lz«l/'^--5«. (6) 

Similar to the treatment of the bands, K98 assumed that the emis- 
sivity of the circumsolar ring was the same as that of the Earth- 
trailing feature, below. We also relax this assumption, and allow 
them to be different. The shape of the expected signal from the 
circumsolar ring is shown in the fourth row of Fig. 4. 
The density of the Earth-traihng feature is given by 

riB (R) = riTB ■ e-(«-«™>'/'^'™-|z«l/'^j«-(«~e™)'/'^2r« . (7) 

We note that for both the circumsolar ring and the Earth-trailing 
feature, there is a typo in the text of K98 - a factor of 2 in the de- 
nominator of the first and third terms in the exponential has been 
added in the text, compared to what is in the code. We follow 
the code. The expected signal from the Earth-trailing feature is 
shown in the fifth row of Fig. 4. 

4.1.4. Integrated Emission 

The total Zodiacal emission is calculated as 

(v) = 6. J t/R ■ (R) Biv,! (R)) , (8) 

where x is the Zodiacal component, v is the frequency, R is a 
location in the Solar System, and B (v, T (R)) is the Planck func- 
tion for the given frequency and temperature at the given loca- 
tion, given by To//?*, with Tq and 6 being parameters, e is the 
emissivity for the given component, which we will be finding 
with our fit. H j is the density for the given component, described 
above. The integral is done along the line of sight, from the lo- 
cation of the satellite to 5.2 AU. 



4.2. Galactic Emission Seen Tlirougli Sidelobes 

Like all real telescopes, Planck has some small sensitivity to 
radiation arriving at the telescope from off-axis, often called 
spillover or 'far sidelobes'. Note that here we are referring to 
light coming from more than five degrees from the instrument's 
nominal pointing direction, which should not be confused with 
'near sidelobes', discussed in the HFI transfer function and 
beams papers (Planck Collaboration VII 2013, see also Tauber 
et al. (2010)). Figure 5 shows the main stray light routes. 



PR Spillover 




Fig. 5. Origin of far sidelobes. The 'SR Spillover' (for 
'Secondary Reflector Spillover'; the lowest set of rays on the 
left of the figure) arrives at the focal plane from outside the sec- 
ondary mirror, directly from the sky. The 'PR Spillover' (for 
'Primary Reflector Spillover'), arrives at the focal plane from 
above the primary mirror and reflects off of the secondary to ar- 
rive at the focal plane. The set of rays between these two contri- 
butions represents the main beam. The 'baffle' contribution, light 
reaching the focal plane after reflecting from the inner sides of 
the baffles, is not shown here. It is often included as part of the 
SR Spillover Adopted from Tauber et al. (2010). 

The secondary reflector (SR) spillover arises from radiation 
that reaches the focal plane without reflecting off of the primary 
reflector. As such, a major component is radiation from the gen- 
eral direction of the telescope boresight, though well outside the 
main beam. The 'baffle' contribution to the SR spillover results 
from radiation which reflects off of the baffles to arrive at the 
focal plane. The primary reflector (PR) spillover arises from ra- 
diation that comes to the satellite from just above the primary 
mirror, reflects off the secondary mirror and arrives at the detec- 
tors. 

At the highest Planck frequencies, the Galactic centre is 
bright enough to be seen through these far sidelobes, even 
if faintly. Since the orientations of these sidelobes change as 
the instrument scans, the survey differences done to detect 
the Zodiacal emission are also sensitive tests of the FSLs. 
Though the Galactic emission mechanism and amplitude is dif- 
ferent, analogous effects are discussed for the LFI in Planck 
Coflaboration IV (2013). 
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Table 1. COBE Zodiacal bands, their associated IRAS band names, and information about the asteroid famihes with which they 
have been associated. 



K98° 


Sc" 


IRAS' 


'p 


Asteroid'' 


'■p 


Type* 


Age" 


a' 


1 


8.78 


7 


9.35 


490 Veritas 


9.26 


C/Ch/... 


8.3 


3.169 


2 


1.99 


a 
P 


1.34 
2.11 


656 Beagle 
832 Karin 


1.34 
2.11 


.../.../ C 
.../.../ S 


<10 
5.8 


3.157 
2.866 


3 


15.0 


J/K 


12.0 


845 Naema 
4562 lannini 


11.96 
12.17 


.../C/... 
.../.../ S 


>20 
<5 


2.644 






M/N 


15.0 


1521 Seinajoki 


15.02 


.../.../... 




2.852 



Notes. The K98 band designation; The 6( parameter of the K98 Zodiacal emission model. This parameter determines roughly at which 
ecliptic latitudes, in degrees, the band appears; The IRAS bands associated with the given COBE band; '''' The modelled proper inclination of 
the given band. The first three are from Grogan et al. (2001), while the last two are from Sykes (1988); The asteroid family associated 
with the given IRAS band. These all come from N03, except that for 832 Karin, which comes from N08; The average proper inclination for 
the associated asteroid family. These all come from N03, except that for 832 Karin, which comes from N08; The spectral type/classifications. 
The first two entries in the triplets correspond to the Tholen and SMASSII classes (Bus & Binzel 2002), while the third corresponds to the SDSS- 
based classification (Carvano et al. 2010). Ellipses are used to indicate that the given classification was not found; The time since the asteroid 
disruption which created the associated asteroidal family, in Myr. These come from N03, except for 832 Karin, which comes from N08; The 
proper semi-major axis of the asteroidal family associated with the band, in AU. These come from N03, except for 832 Karin, which comes 
from N08. 



To this point, resource constraints have limited this study to a 
single far sidelobe calculation for all detectors. We use a GRASP 
calculation of the far sidelobes for the 353-1 horn (see Fig. 9 of 
Planck Collaboration VII 2013) and do not attempt to correct 
for differences in frequency or location for other horns. While 
this is not optimal, and will be improved in later releases, we 
note that the primary, large-scale features of the far sidelobes are 
defined by the telescope, rather that the horns or their placement, 
so first-order effects should be captured. Some of the limitations 
imposed by this are discussed in Sect. 6.1. 

To make templates of what we might see from the Galactic 
centre through the far sidelobes, we use the Planck simulation 
software described in Reinecke et al. (2006), but use GRASP 
calculations of the the various far sidelobe components instead 
of main beams. As inputs to the simulations, we use the actual 
Planck maps at the appropriate frequency. This allows us to ac- 
count for differences in the brightness of the Galaxy as a function 
of frequency with minimal effort - the Galactic templates should 
already be calibrated in the correct units at the given frequency. 
The far sidelobe templates are made at the timeline level and run 
through the relevant parts of the pipeline software. In particular, 
the offset removal, or 'destriping' must be done on the templates 
before fitting in order to get reasonable fit results. The result- 
ing templates made using these FSL calculations with a Planck 
857 GHz sky as input are shown in the bottom three rows of 
Fig. 6. 

One factor for which we do not account with these tem- 
plates is the difference in spillover between the different frequen- 
cies (Lamarre et al. 2010; Tauber et al. 2010). Since we illumi- 
nate more of the telescope at lower frequencies than at higher 
frequencies, if our templates simply scale with the spillover, we 
will expect different fit values for our templates at different fre- 
quencies. 

To illustrate the relative contributions of these various tem- 
plates, in Fig. 7 we show a series of maps. Maps in each row are 
similar to those in the previous row, except that one more tem- 
plate or group of templates has been added to form the new row. 
For all rows, the first column corresponds to data from the first 
Planck survey, the second column corresponds to data from the 
second Planck survey, and the third column is the difference of 
the second column minus the first. The first row shows the sum 
of all the far sidelobes; the second row shows the result when we 
add dust band 1 to the far sidelobes - note that the scales change 



from the first to the second row; the third row shows the sum of 
the far sidelobes and the first two bands; the fourth row shows 
the sum of the far sidelobes and all the dust bands; the fifth row 
shows what's in the fourth row, plus the circumsolar ring and 
Earth-trailing feature; finally, the bottom row shows the sum of 
the far sidelobes and all Zodiacal templates (and the scales have 
again changed). As this is simply illustrative, we have assumed 
unit emissivities for the Zodiacal components, and multiplied the 
far sidelobe components by a factor of 15 (which we will see in 
Sect. 5 is representative of the most extreme, multi-moded case). 
With these caveats in mind, the survey difference of the sums of 
all the components, the lower right image, can be compared to 
the bottom image in Fig. 3. 

5. Spectrum 

We fit the data shown in the bottom panel of Fig. 3, as well as 
the analogous data at other frequencies and for the second year 
of observations, to a constant plus combinations of the templates 
shown in Sect. 4, Figs. 4 and 6. 

For each fit, we assume that the survey difference map at sky 
pixel p, called Dp, can be modelled as 

Dp^^ (e2,/ ■ T2,t,p - £i,t ■ + constant, (9) 

t 

where e[i|2],/ is the emissivity fit for template t during survey 1 
or 2 at the given frequency and r[i|2j,/,p is the value of the f* 
template at pixel p for survey 1 or 2, calculated as described in 
Sect. 4. For example, for our 'basic' fit, we will have 19 tem- 
plates - one each for the diffuse cloud, circumsolar ring. Earth- 
trailing feature and each of the three dust bands, as well as one 
each for the Galactic far sidelobes. All of these are repeated 
twice; once for each survey in a yearly difference map. Finally 
we also fit to an overall constant, to which Planck is not sensi- 
tive. We then minimize 

x'-T,{^^-Dpf■ (10) 

p 

Separating the templates into surveys has the disadvantage of 
increasing the number of parameters in our fits. While we do not 
expect either the emissivity of the Zodiacal emission or the far 
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Fig. 6. Templates of the Dipole and the Galaxy seen through our far sidelobes for sui-veys 1 (left) and 2 (center), and for the difference 
in these two (right). Row 1: Dipole seen through the direct SR contribution; Row 2: Dipole seen through the PR contribution; Row 
3: Dipole seen through the baffle SR contribution; Row 4: Galaxy seen through the direct SR contribution; Row 5: Galaxy seen 
through the PR contribution; Row 6: Galaxy seen through the baffle SR contribution; These last three templates are made using 
857 GHz data passed through the far sidelobe calculation described in Sect. 4.2. The scales are different for the top and bottom three 
rows, and for the first two and the last column. 
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Fig. 7. A sequence of maps, each building upon that above, designed to show the relative contributions of the various templates in 
Figs. 4 and 6. Row 1: The sum of all templates in Fig. 6; Row 2: Row 1 plus dust band 1; Row 3: Row 2 plus dust band 2; Row 
4: Row 3 plus dust band 3; Row 5: Row 4 plus the circumsolar ring and Earth-trailing feature; Row 6: The far sidelobes and all 
Zodiacal components. Left: Survey 1; Center: Survey 2; Right: Survey 2 minus Survey 1. Note that the amplitudes used are only 
approximate, the figure being for illustrative purposes only. The scales change between the first and second rows, and between the 
fifth and sixth rows, as well as between the second and third column. 
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sidelobe calculations to change from survey-to-survey, we sep- 
arate them in this way for two reasons. The first is simply as a 
basic reality check - if we see significant differences between fits 
to two different surveys, we should be sceptical. Beam asymme- 
tries or transfer function effects, for example, might cause dif- 
ferences from survey-to-survey, as might imperfections in the 
model itself. 

The second reason is to calculate the error bars. As just 
noted, we are often as concerned by systemic effects as much 
as by "random" noise. By separating the data by survey, we may 
calculate error bars using the standard eiTor of the successive 
measurements as a proxy for the uncertainties, rather than prop- 
agating white noise estimates. This should allow us a more con- 
servative estimate of our uncertainties, which accounts for model 
deficiencies or low levels of systematics which change by survey 
(the aforementioned beam asymmetries and transfer functions, 
for example). 

A note about weighting: the bulk of the Zodiacal emission 
is, of course, in the ecliptic plane. Planck, on the other hand, has 
more statistical weight, on a pixel-by-pixel basis, at the ecliptic 
poles. We therefore use uniform weights over pixels, rather than 
statistical weights, since this would down-weight specifically the 
regions with our signal. 

As mentioned above, we fit each of the two years' survey dif- 
ference maps to a cloud, circumsolar ring. Earth-trailing feature, 
three bands and three far sidelobe templates, plus a constant. 
The results for the four 857 GHz horns are shown in Fig. 8. 
Averaging over horns and surveys at all six HFI frequencies 
yields Fig. 9. Numerical values are given in tables 2, 3 and 4. 



Table 2. Emissivities of the diffuse cloud, circumsolar ring, and 
Earth-trailing feature from fit result averages. 



Table 4. Fit coefficients for the Galaxy seen through the far side- 
lobes. 



V (GHz) 


Cloud 


Rin 


g 


Feature 


857 


0.301 


± 0.008 


0.578 ± 


0.359 


0.423 ± 


0.114 


545 


0.223 


± 0.007 


0.591 ± 


0.203 


-0.182 ± 


0.061 


353 


0.168 


± 0.005 


-0.211 ± 


0.085 


0.676 ± 


0.149 


217 


0.031 


± 0.004 


-0.185 + 


0.143 


0.243 ± 


0.139 


143 


-0.014 


± 0.010 


-0.252 ± 


0.314 


-0.002 ± 


0.180 


100 


0.003 


± 0.022 


0.163 ± 


0.784 


0.252 ± 


0.455 



Table 3. Emissivities of the three dust bands from fit result aver- 
ages. 



V (GHz) 


Band 1 


Band 2 


Band 3 


857 


1.777 ± 


0.066 


0.716 


± 0.049 


2.870 ± 


0.137 


545 


2.235 ± 


0.059 


0.718 


± 0.041 


3.193 + 


0.097 


353 


2.035 ± 


0.053 


0.436 


± 0.041 


2.400 ± 


0.100 


217 


2.024 ± 


0.072 


0.338 


± 0.047 


2.507 ± 


0.109 


143 


1.463 ± 


0.103 


0.530 


+ 0.073 


1.794 ± 


0.184 


100 


1.129 ± 


0.154 


0.674 


+ 0.197 


1.106 + 


0.413 



6. Discussion 

As discussed in the HFI processing paper (Planck Collaboration 
VI 2013) and in the Planck explanatory supplement Planck 
Collaboration ES (2013), these fits are used to create the im- 
plied Zodiacal and sidelobe emission in each HFI observation, 
which can then be removed and the maps are recreated. The sur- 
vey 2 minus survey 1 difference maps for the 857-1 horn both 



V 

(GHz) 


PR Spill. 


Pred." 


SR Spill 
Direct 


over 
BafHe 


Pred.* 


100 


-25.8 ± 5.7 


7 


26.3 + 15.7 


56.7 + 5.7 


10 


143 


-9.1+4.1 


6 


13.0 + 4.8 


23.8 + 5.4 


10 


217 


0.6+ 1.1 


5 


6.3 + 2.1 


6.5 + 1.3 


6 


353 


-1.2 + 0.5 


1 


-4.3 + 2.1 


3.6 + 0.7 


1 


545 
857 


7.7 + 1.7 
17.1 + 3.4 


15 
1.5 


8.8 + 3.1 
23.9 + 4.2 


7.9+ 1.0 
16.7 + 3.1 


1 

0.005 



Notes. The (unit-less) value we would expect for the fit to the 
primary spillover sidelobe contribution. It is the ratio of the primary 
spillover at the given frequency over the spillover at 353 GHz, as cal- 
culated in table 2 of Tauber et al. (2010). If all our data and predictions 
were perfect, this value would match the corresponding fit value in the 
Col. labelled 'PR Spill.'; The (unit-less) value we would expect for 
the fit to the secondary spillover sidelobe contribution. It is the ratio 
of the secondary spillover at the given frequency over the spillover at 
353 GHz, as calculated in table 2 of Tauber et al. (2010). If all our data 
and predictions were perfect, this value would match the corresponding 
fit values in the Cols, labelled 'SR Spill.'. 



before and after Zodiacal and far sidelobe removal, are shown in 
Fig. 10. In this section we discuss the fit impUcations further. 

6.1. Far Sidelobes 

In addition to fit values obtained for the Galaxy seen in each 
component of the far sidelobes, the Cols, labelled 'Pred.' in ta- 
ble 4 show the expected values of the spillover, normalized to 
that of the 353 GHz channel, from Tauber et al. (2010). These are 
the ratios of the expected spillover in each frequency, compared 
to that at 353 GHz, the frequency for which the sidelobe calcu- 
lations were done. Since the fit values account for the changes 
in Galactic emission with frequency, if our predictions and data 
were perfect, the fit values would match those of the predictions. 

The FSL signature is clearly visible at 857 GHz in the bot- 
tom panel of Fig. 3, and this is reinforced in table 4. As the 857 
and 545 GHz channels are multi-moded, the differences are not 
that surprising; it is difficult to do the calculations necessary for 
the prediction. In addition, the specifications for the horn fabri- 
cation were quite demanding, and small variations, though still 
within the mechanical tolerances, could give large variations in 
the amount of spillover. 

For the lower frequency, single-moded channels, however, 
the situation is different. There is no clear detection of PR 
spillover. While the significant negative values may indicate 
some low-level, large-scale systematic, there seems to be noth- 
ing with the distinctive signature of primary spillover at frequen- 
cies between 100 and 353 GHz. 

For the direct contribution of the secondary SR spillover, 
the situation is similar at 353 GHz, but at 217 and 143 GHz we 
are finding a 3cr detection at about the level expected, while 
at 100 GHz the value is about 2.5 times higher than expected, 
though the signal-to-noise on the detection is less than 2cr. The 
baffle contribution to the SR spillover seems to be in accord with 
expectations at 353 and 217 GHz, and higher than what is pre- 
dicted at 100 GHz. 

The values for the PR spillover, which is the most distinctive 
of the far sidelobe patterns, and presumably the easiest to dis- 
entangle from other effects, would indicate that the PR spillover 
values in table 2 of Tauber et al. (2010) may be slightly overesti- 
mated. The values for the direct contribution of the SR spillover 
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Fig. 8. Results of fits at 857 GHz. Two sets of fits are shown in each panel; one to the first year of Planck observations, in red, and 
one to the second year of Planck data, in blue. For each of these fits, the emissivities found for surveys 1 and 3 are shown with 
circles, while surveys 2 and 4 are shown with squares. In the absence of time variability of the Zodiacal emission, we would expect 
little difference between coiTesponding red and blue symbols. Agreement or disagreement between squares and circles give some 
indication of the systematics in the data and of the of the coiTectness of the templates. The average of all four measurements for each 
horn is shown as the black square, and the average of all red and blue circles and squares for the entire frequency, and their standard 
errors, are represented with the horizontal gray band. Where appropriate, a dotted line is used to mark the zero level. Similar plots 
for all HFI frequencies can be found in the Planck Explanatory Supplement (Planck Collaboration ES 2013). 



roughly confirm the far sidelobe calculations. The baffle contri- 
bution to the SR spillover seems a bit high. We take the ensem- 
ble of these numbers as rough confirmation that our beam calcu- 
lations are not drastically incorrect, but do not use the specific 
numbers in either Planck Collaboration Vll (2013) or Planck 
Collaboration VIII (2013). Similar conclusions are drawn for the 
LFI in Planck Collaboration IV (2013). 

We have also included a template of the Dipole as seen 
through the far sidelobes in the fits in some fits to check to see if 
they are detected. As expected, they are not. The results quoted 
above are from fits which do not include these Dipole templates. 



tional to the frequency-squared at longer wavelengths, for com- 
parison with Fig. 2 of Fixsen & Dwek (2002), who used FIRAS 
data to investigate the far-infrared/sub-millimetre behaviour of 
the Zodiacal cloud. These results are consistent with their con- 
clusions^. 

As the diffuse cloud is so much brighter in the mid-infrared 
than in the CMB, its relatively low level at Planck wavelengths 
has been exploited in Planck Collaboration IX (2013) to set in- 
flight limits on any possible out-of-band leaks in the instrument's 
spectral transmission. 



6.2. Diffuse cloud 

Figure 9 shows the emissivity of the diffuse cloud falling off 
with increasing wavelength, as would be expected for particles 
with characteristic sizes of order 30 fim. The dashed line in this 
figure shows a flat emissivity to 150 fim, with values propor- 



Care must be taken with direct comparisons between the emissiv- 
ities quoted in Fixsen & Dwek (2002), who quote emissivities relative 
to a 245 K cloud, and those in this work, which follow K98 and assume 
a cloud of temperature 286 K at 1 AU from the Sun. 
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Fig. 9. Emissivities of various components of the K98 Zodiacal emission model obtained from both Planck/HFl and COBE/DIRBE. 
The COBE/DIRBE points, at wavelengths of 250 microns and less, are shown to the left of the plot over the grey background. The 
Planck/HFl points are to the right of the plot, at wavelengths greater than 250 micron. In both cases, the diffuse cloud is shown as 
blue squares. For DIRBE, the dust bands were assumed to all have the same emissivities, shown as the red, left-pointing triangles. 
For HFI, the bands were allowed to have different emissivities, which are shown as red, up-pointing triangles for dust band 1, as 
pink, right-pointing triangles for dust band 2, and as orange, down-pointing triangles for dust band 3. Similarly, K98 assumed that 
the circumsolar ring and Earth-trailing feature had the same emissivities, shown as the green circles. For Planck we allow them to 
be different, shown with the green and yellow hexagons, respectively. The Planck values were obtained by fitting for an amplitude 
to each of these, as well as Galactic sidelobes. All other parameters in the model were assumed to have their K98 values. Each point 
is the average of the coiTesponding values obtained for all individual horns and surveys with the given frequency, including the first 
two years of HFI data - the gray shaded regions in each panel of Fig. 8. The error bars shown are from the standard eiTors of these 
different measures. The numerical values are given in tables 2 and 3. Note that a few cloud, circumsolar ring, and Earth-trailing 
feature values are negative, and so do not appear in this log-log plot. In such cases, the upper limit will appear as a short horizontal 
line. The dotted line indicates an emissivity of unity at all wavelengths, and the dashed line indicates an emissivity which is unity at 
wavelengths below 150 //m and which is inversely proportional to wavelength-squared at longer wavelengths. 



6.3. Circumsolar ring and Eartli-trailing feature 

We draw no conclusions about the circumsolar ring or the Earth- 
trailing feature. The fit values obtained for their emissivities are 
inconsistent from frequency-to-frequency, their values often be- 
ing negative. This remains true for Planck data when the fits are 
done requiring that the two components have the same emissiv- 
ities, as was done in K98. 

Inspection of the panels labelled 'Ring' and 'Feat.' (i.e., the 
middle- and lower-left panels) of Fig. 8 show systematic dif- 
ferences between results from even- and odd-numbered surveys 
(that is, the circles and squares seem to be systematically dif- 



ferent, regardless of whether they are blue or red). Planck's ob- 
serving pattern was different for odd and even numbered sur- 
veys, but very similar for either even numbered surveys alone, or 
for odd numbered surveys alone. Random noise is not an issue, 
since the measurements with similar observations are repeat- 
able, so this is either an indication that the circumsolar ring and 
Earth-trailing model templates themselves need improvement at 
Planck wavelengths, or that systematic errors are affecting these 
specific components. 
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Fig. 10. Survey 2 minus survey 1 difference maps before (left) and after (center) Zodiacal emission removal, as well as the total 
(i.e., not differenced) Zodiacal emission removed from the nominal mission HFI maps (right). The rows, from top to bottom, are for 
857, 545, 353, 217, 143, and 100 GHz. The top two rows, for 857 and 545 GHz, are in units MJy/sr, while the other frequencies are 
in units of yuKcMB- Note that while these images are presented together for compactness, they are different, as some are differences 
and some or not. 
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We note that the conclusions presented elsewhere in this 
work remain essentially the same whether or not we include the 
circumsolar ring and Earth-traihng feature. 

6.4. Bands 

An interesting feature of Fig. 9, and the primary result of this 
work, is the difference in the emissivities of the bands compared 
to that of the diffuse cloud, pointing towards differences in the 
particles in these different components. While there may be hints 
of this in the longest wavelength DIRBE data, the effect becomes 
clear at Planck wavelengths. This is not necessarily unexpected. 
The composition of the diffuse cloud is still disputed, but is often 
claimed to be both asteroidal and cometary (see, for example, 
Kortenkamp & Dermott 1998; Nesvorny et al. 2010; Tsumura 
et al. 2010). Since the bands, on the other hand, are understood 
to be asteroidal debris only (Sykes & Greenberg 1986), the dif- 
ference may simply be a reflection of these different origins. 

The fact that the fitted emissivities of bands 1 and 3 rise 
above unity is perplexing. At first glance, one might imagine 
some new, cold component in the cloud which would cause an 
excess that might be interpreted as an excess in emissivity in 
some other component. However, to peak around 545 GHz, this 
component would have to have a temperature of the order of 
lOK, and therefore be much more distant than most of the dust 
usually associated with the Zodiacal cloud. While enticing, it is 
difficult to understand how such a component could survive the 
differencing process used in this analysis, which reduces signals 
from distant sources more than those nearby, or how such a com- 
ponent could mimic an excess in two dust bands above from the 
ecliptic plane, but not do the same in the other Zodiacal cloud 
components. 

One might worry that covariance between between the var- 
ious components might be causing problems in the fitting pro- 
cedure. To check this, we have repeated the fit including and 
omitting various combinations of the the circumsolar ring and 
Earth-trailing feature, or both, and assuming their emissivities 
were independent or equal. In no case did the difference between 
the diffuse cloud and the dust bands disappear. 

The excess may ultimately be explained by degeneracies 
in the model for the density of the bands. As presented in 
Sect. 4.1.2, the normalisation of the density of particles is com- 
pletely degenerate with the emissivity for each band. In addition, 
the emission is also roughly proportional to the temperature nor- 
malisation, because we are observing in the Rayleigh-Jeans tail 
of the Zodiacal emission. While any overall change in the tem- 
perature of the IPD particles would scale all components of the 
Zodiacal emission, because temperature is nearly inversely pro- 
portional to distance from the Sun, the location of the bands are 
important. While the excess over 1 is too large to be explained by 
errors in distance and thus temperatures alone, one might appeal 
to a change in a combination of distance, particle density nor- 
malisation and emissivity of these bands in particular to arrive at 
mutually consistent results for both Planck and DIRBE. As this 
will involve simultaneous work with both Planck and DIRBE 
data, it is beyond the scope of this paper. 

Bands 1 and 3 also seem to show different behaviour than 
band 2. Since bands 1 and 3 are both at high ecliptic latitude, 
while band 2 is not, one might again worry that one of the other 
templates to which we are fitting might have significant overlap 
with a subset of the bands, which in turn causes an apparent 
difference in emissivities. To check this, we have repeated the fits 
with and without various combinations of the cloud, circumsolar 
ring and Earth-trailing feature, as well as the far sidelobes. In all 



cases, bands 1 and 3 are always significantly different than dust 
band 2. When the diffuse cloud itself is omitted from the fit, 
the emissivity of dust band 2 goes up, but is still distinctively 
different from that of dust bands 1 and 3. 

As dust band 2 is a combination of the IRAS a and fi bands, 
one may also worry that one of these two is more important 
for the shorter IRAS and COBE wavelengths, but that the other 
might be more important for the longer Planck wavelengths. We 
note that its emission is dominated by contributions from the 
Karin/Koronis family (see Nesvorny et al. 2008, Fig. 1), but have 
therefore confirmed specifically that varying the parameter of 
the second band between values appropriate for either a or /3 
does not remove this difference (see table 1). 

If the age of band 2 were significantly different from those 
of bands 1 and 3, we might argue that Poynting-Robertson drag 
had depleted some of the bands of more smaller particles than 
the others (Wyatt et al. 2011). N03 and N08, however, have es- 
timated the ages of most of the asteroid families which might be 
associated with the bands (reproduced in table 1), and the age of 
any of the associations with band 2 is between those of any of 
the possible associations with bands 1 or 3. These same figures 
tend to rule out modifications of the material properties due to 
photo-processing or solar wind exposure for differing periods. 
Band 2 also seems to be roughly the same distance from the Sun 
as the other two bands, so it is difficult to appeal to differences 
in environment as the cause. 

We speculate on the following to explain any differences: 
Veritas, the asteroid family proposed to be associated with dust 
band 1, is classified as carbonaceous (Bus & Binzel 2002). 
As noted above, the IRAS /3 band, associated with the Karin 
family of asteroids, seems to dominate the emission from dust 
band 2. Karin and its larger sibling, Koronis, are classified as 
siliceous, or stony, objects (Bus & Binzel 2002; Carvano et al. 
2010). While the dust band 3 has a number of asteroid fami- 
lies which may be contributing to it (see table 1), we propose 
that the emission is dominated by carbonaceous-based asteroid 
families (three quarters of the asteroids in the Solar System are 
carbonaceous), and that the difference in emissivities between 
dust band 2 and dust bands 1 and 3 arises from this difference in 
composition. The differing emissivities may be either due to this 
intrinsic composition difference, or the size-frequency distribu- 
tion of particles that results from different kinds of asteroids col- 
liding (Grogan et al. 2001, for example). This explanation would 
not be valid if it were to turn out that dust band 3 was dominated 
by dust associated with the lannini asteroid, for example, as it is 
siliceous. 

6.5. Implications for the CMB 

The right-hand column of Fig. 10 shows the Zodiacal emission 
implied by the fits, created by subtracting the maps made apply- 
ing the Zodiacal emission correction from those that were made 
without it applied. One can see here the difference in the rela- 
tive amplitudes of the bands versus the diffuse cloud, the bands 
being relatively more important at low frequencies than at high. 
Figure 11 shows the power spectra of the Zodiacal correction 
maps, all in units of {jjKcmb)^- Here, the cloud is seen at multi- 
poles of less than about 10, while the bands and other structures 
are see in higher multipoles. 

At 143 GHz, The signal reaches a few //Kcmb in the 
map, while the power spectrum has values of the order of 
one (yuKcMfi)^- The absence of power in the odd multipoles is 
a reflection of the north-south symmetry of the signal. These 
can be compared with the CMB temperature spectrum, which 
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Fig. 11. Power spectra of the Zodiacal correction maps shown in 
the right-hand column of Fig. 10. Black plus signs, grey crosses, 
cyan triangles, orange diamonds, blue circles, and red squares 
represent, respectively, 857, 545, 353, 217, 143, and 100 GHz. 
For the CMB channels, even multipoles are shown with filled 
symbols, while odd multipoles are shown with empty symbols. 
This "even-odd" pattern is a reflection of the symmetry around 
the ecliptic plane - odd multipoles are almost absent, as they 
would indicate structure in the maps that were anti-symmetric 
about the ecliptic plane. A standard ACDM CMB temperature 
spectrum is shown as the grey line roughly half-way down the 
plot, orders of magnitude above the Zodiacal spectrum in the 
Planck CMB channels. 

is shown as the black line about halfway between the top and 
the bottom of the plot^. The Zodiacal emission correction spec- 
tra are orders of magnitude smaller than the CMB spectrum. The 
Zodiacal emission does not compromise Planck's cosmological 
results. 



7. Conclusion 

Zodiacal emission has long been an important foreground for 
searches for the extra-galactic background at infrared wave- 
lengths. With the ever-increasing sensitivity of CMB experi- 
ments, it will soon become important to account for in sub- 
millimetre and CMB analyses as well. 

The K98 model does fairly well in modelling the diffuse 
Zodiacal cloud emission at Planck wavelengths, as long as ap- 
propriate emissivities are assumed. It does less well, however, 
in modelling the other features. Because they appear to be more 
emissive at these frequencies than the cloud, the bands contribute 
more to the Zodiacal emission relative to the diffuse cloud at 
CMB frequencies (i.e., near 143 GHz, or 2.1mm). The 2013 
Planck release includes both maps which have had Zodiacal 
emission removed, and maps which have not had Zodiacal emis- 
sion removed. We note that Planck Collaboration Xll (2013) 
found better results using maps that had not had the Zodiacal 
emission removed to make estimations of the CMB and other 
astrophysical components in the HFl maps. The component sep- 
aration methods used there naturally correct for a large amount 
of Zodiacal emission, as it is spectrally similar to Galactic dust 
emission in the Planck CMB channels. In addition, one can see 
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in Fig. 9 and table 2 that the diffuse cloud is not actually de- 
tected at 143 GHz, the frequency which carries the most weight 
in Planck CMB analyses. The Zodiacal emission removal, of 
course, depends delicately on how this non-detection is treated. 
This will be addressed in future Planck releases. 

Improvements in modelling of the circumsolar ring. Earth- 
trailing feature, and the dust bands, as well as inclusion of fainter 
and partial bands (e.g.. Espy et al. 2009) should be done to make 
truly "clean" CMB maps. As these bands are believed to be 
the products of asteroid collisions, further study of these bands 
at these wavelengths may also inform us on the nature of the 
intermediate-sized particles created during the destruction of the 
associated asteroids. We may hope to learn not only more about 
the size distribution, but also the differences between, for exam- 
ple, the results of a collisions involving siliceous and carbona- 
ceous asteroids. 

Just as material from asteroid collisions contributes to 
Zodiacal emission, material shed from comets must also con- 
tribute. As part of the HFI data reduction process, we mask out 
solar system objects which would cause "noise" in the final sky 
maps. We searched for comets as part of this process, but found 
only one (Christensen). For completeness, appendix A presents 
this information, as it is not possible to extract from the data in 
the 2013 Planck map release. We have not yet detected extended 
tails of comets. 

One of the primary goals of the next stage of analyses, once 
this "nearby" IPD has been completely removed, will be to 
search for dust associated with the Kuiper Belt, or set limits upon 
it. 

The full-mission Planck data release will include polarisa- 
tion information. While it is not expected, limits will be put on 
possible contamination of the polarisation of the CMB by po- 
larised Zodiacal emission, as has been done here for tempera- 
ture. 

Work is now underway to address all these points for the next 
Planck data release. While the signal is quite small - at CMB 
wavelengths the signal we are discussing is orders of magnitude 
smaller than the primary CMB anisotropics, it is detectable and 
should be subtracted from the data. There will be improvements 
in dust modelling, improvements in satellite modelling, and ad- 
ditions to address polarisation. The ultimate goal will be simulta- 
neous analyses with IRAS, COBE, and Akari and other to under- 
stand the large-scale Zodiacal emission from the near-infrared to 
the microwave. 
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Appendix A: Solar System Objects 

HFI detects a number of moving solar system objects; mainly 
asteroids but also one comet. Most of these are masked at an 
early stage of the analysis so as not to affect the maps of the sky. 
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This, however, means that they cannot be easily extracted from 
the delivered products. We therefore present them here. 

Their location is computed using the JPL Horizons^ 
(Giorgini et al. 1996) system programmed with Planck's, orbit. 
Table A.l shows the epochs of Planck's, observation of these ob- 
jects and tables A. 2 and A. 3 show the distance from the Sun and 
Planck at the time of observation, as well as the position on the 
sky. 

During the standard HFl timeline processing, these objects 
are flagged and not included in the standard HFI maps (Planck 
Collaboration VI 2013), which makes these maps an excellent 
tool for removing the background of these moving objects. 

5Astraea, Christensen, and 128Nemesis were not flagged in 
the maps, so their fluxes reported here have been adjusted for 
the fact that a fraction of their emission would have still been in 
the maps used to remove the background before estimating their 
fluxes. 

We select time-ordered data samples within 0?5 of the 
source and recalibrate into MJy/sr (using IRAS-conventions). 
We project the pointing of Planck into coordinates relative to 
the predicted position of the moving object. A synthetic back- 
ground timeline is estimated by re-sampling the HEALPix- 
gridded (Gorski et al. 2005) Planck maps using cubic spline in- 
terpolation. 

We fit the main beam template (Planck Collaboration Vll 
2013) for each bolometer to the time-ordered data. There are 
seven free parameters in the fit: xq and jq, corresponding to the 
centroid of the object, a rotation angle i/f, an amplitude A, and 
three parameters describing a linear slope in x and y of any resid- 
ual background. The amplitude times the solid angle of the beam 
model gives the flux density. Figure A. 1 shows an example. 

We find negligible difference between a fit assuming a 
Gaussian template for the beam instead of the PSF, as reported 
in Planck Collaboration Vll (2013), so here we report the PSF 
fit flux densities in table A. 4. We also tried aperture photome- 
try, but the results were noisier and inconsistent from season-to- 
season, which might be expected, as the Ecliptic plane is some- 
what under-sampled when using only a single season of data. 

A.1. Notes 

The residual map variance can vary quite a bit due to the very 
different backgrounds. During season one, for example, 1 Ceres 
was in a region of very high foregrounds, and so was difficult 
to detect at 545 and 857 GHz, the bands most susceptible to 
foregrounds. 

A.2. Basic Behaviour 

Figure A. 2 demonstrates some of the basics of the asteroid de- 
tections. For each measurement, we denote the measurement as 
/, the distance between Planck and the object as d, and the dis- 
tance between the Sun and the object as s. These quantities for 
the first measurement of each asteroid will be denoted fit, do, 
and So- In the top panel, we show asteroids that were detected in 
multiple surveys. Assuming the temperature of the object at any 
time is inversely proportional to V^, we would expect 

fV~s^^, (A.l) 
which is roughly seen in the data. 
^ http : //ssd . jpl ■ nasa . gov/?hor izons 



In the bottom panel, we show the 545-to-857 GHz spectral 
indices, as well as the 545-to-353 GHz spectral indices for Ceres 
and Vesta. 
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Fig. A.l. An example of background removal and the residual after PSF fitting for the second observation season of ICeres at 
857 GHz. Background removal and fitting are done in the time domain, bolometer by bolometer; this stacked map of all the 857 GHz 
data is created only for visualization. Left: Stacked data before background removal; Center: Stacked data after background removal; 
Right: Stacked data after background and source removal. 



Table A.l. Observation seasons of Solar System objects detected by Planck. 



Object Season 1 Season 2 Season 3 Season 4 Season 5 



ICeres 


24/03/10- 


-25/03/10 


10/09/10- 


-11/09/10 


03/07/11- 


-04/07/11 


03/12/11- 


-04/12/11 






2Pallas 


23/02/10- 


-25/02/10 


30/07/10- 


-31/07/10 


02/05/11- 


-03/05/1 1 


13/10/11- 


-15/10/11 






3Juno 


10/12/09- 


-12/12/09 


30/12/10- 


-01/01/11 


27/05/11- 


-28/05/1 1 










4Vesta 


24/11/09- 


-26/11/09 


29/04/10- 


-01/05/10 


28/04/11- 


-29/04/11 


26/10/11- 


-28/10/11 






5Astraea 


05/10/09- 


-06/10/09 


30/07/10- 


-31/07/10 


02/01/11- 


-03/01/11 


27/12/11- 


-28/12/11 






6Hebe 


27/06/10- 


-30/06/10 


13/12/10- 


-15/12/10 


21/12/11- 


-22/12/11 










SFlora 


12/06/10- 


-15/06/10 


30/1 1/10- 


-02/12/10 


06/01/12- 


-07/01/12 










9Metis 


30/01/10- 


-31/01/10 


11/07/10- 


-13/07/10 


23/04/11- 


-25/04/11 


10/10/11- 


-11/10/11 






lOHygiea 


12/11/09- 


-14/11/09 


18/04/10- 


-19/04/10 


23/02/11- 


-23/02/11 


04/08/11- 


-05/08/11 






1 IParthenope 


01/10/09- 


-02/10/09 


10/03/10- 


-11/03/10 


23/01/11- 


-25/01/11 


03/07/11- 


-05/07/11 






12 Victoria 


20/02/10- 


-21/02/10 


15/08/10- 


-16/08/10 


09/09/11- 


-11/09/11 


21/09/11- 


-21/09/11 






1 3Egeria 


21/03/10- 


-22/03/10 


09/09/10- 


-10/09/10 


17/07/11- 


-18/07/11 


24/12/11- 


-25/12/11 






ISEunomia 


29/03/10- 


-30/03/10 


15/09/10- 


-16/09/10 


27/08/11- 


-28/08/11 










16Psyche 


18/10/09- 


-19/10/09 


16/09/10- 


-17/09/10 


01/03/11- 


-01/03/1 1 


23/12/11- 


-24/12/11 






17Thetis 


03/09/09- 


-04/09/09 


12/02/10- 


-13/02/10 


30/11/10- 


-01/12/10 


02/05/11- 


-03/05/11 






1 SMelpomene 


21/01/10- 


-23/01/10 


06/01/11- 


-07/01/11 


02/06/11- 


-04/06/11 










19Fortuna 


14/09/09- 


-16/09/09 


02/03/10- 


-03/03/10 


26/01/11- 


-27/01/11 


02/07/11- 


-04/07/11 






20Massalia 


13/12/09- 


-15/12/09 


28/12/10- 


-30/12/10 


03/06/11- 


-05/06/11 










29Amphi trite 


02/04/10- 


-03/04/10 


21/09/10- 


-22/09/10 


10/08/11- 


-11/08/11 










41 Daphne 


27/11/09- 


-28/11/09 


28/08/10- 


-28/08/10 


03/02/11- 


-04/02/11 


05/11/11- 


-07/11/11 






45Eugenia 


19/04/10- 


-20/04/10 


06/10/10- 


-07/10/10 


06/08/11- 


-07/08/11 


11/01/12- 


-12/01/12 






52Europa 


23/09/09- 


-24/09/09 


05/03/10- 


-06/03/10 


14/01/11- 


-15/01/11 


21/06/11- 


-23/06/1 1 






88Thisbe 


08/11/09- 


-09/11/09 


25/09/10- 


-26/09/10 


03/03/11- 


-03/03/1 1 


06/12/11- 


-06/12/11 


16/12/11- 


-16/12/11 


128Nemesis 


24/09/09- 


-25/09/09 


05/03/10- 


-06/03/10 


05/01/11- 


-07/01/11 


06/06/11- 


-08/06/11 






324Bamberga 


03/10/09- 


-04/10/09 


12/03/10- 


-14/03/10 


02/01/11- 


-04/01/11 


25/05/11- 


-26/05/11 






SUDavida 


19/02/10- 


-20/02/10 


31/07/10- 


-01/08/10 


01/04/11- 


-02/04/11 


07/09/11- 


-08/09/11 






704Interamnia 


18/02/10- 


-19/02/10 


03/08/10- 


-04/08/10 


16/04/11- 


-17/04/1 1 


01/10/11- 


-02/10/11 






Christensen 


29/09/09- 


-30/09/09 


18/04/10- 


-19/04/10 


11/09/10- 


-11/09/10 


23/03/11- 


-24/03/11 


19/08/11- 


-20/08/11 
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Table A.2. Parameters for observation seasons of Solar System objects detected by Planck. 



Object 


Season 


Solar Range 


Planck Range 


Ecliptic Longitude 


Ecliptic Latitude 






[AU] 


[AU] 


[deg] 


[deg] 


J- C& 


I 


2.774 


2.536 


271.1 


2.169 


1 Ceres 


2 


9 800 


2.624 


963 4 


-4 669 


J- V^Cl Co 


3 


2.982 


2.595 


359.3 


-11.99 


1 Ceres 


4 


9 038 




348 9 


-11.24 


z.jr ail do 


I 


2.632 


2.233 


231.9 


26.18 




2 


2.994 


2.807 


218 


35.55 




3 


3.377 


3.244 


313.5 


32.56 




4 


^ 40^ 


J, i-j 1 


907 


24.6 


3 Juno 


I 


2.051 


1.667 


356.8 


-9.549 


3 Juno 


2 


2.588 


2.21 


178.5 


-3.659 


3 Juno 


3 


2.959 


2.615 


165.4 


3.951 


*+ VColct 


J 


2.472 


9 936 


151 6 


3 140 


4Vest3 


2 


2.325 


1.872 


143.3 


7.996 


4 Vests, 


1. 

J 


2.174 


1 0'S9 


310 1 


916^ 


4Vesta 


4 


2.318 


1.966 


311.1 


-6.731 


S A ctrapji 




3.065 


2.7 


294.2 


0.9236 


S A ctTapa 


2 


2.877 


2.461 


24.35 


-4.22 


S A ctfapn 


3 


2.601 


2.376 


13.72 


-5.692 


S A Qtfapa 


4 


2.085 


1.634 


173.1 


0.4027 


UXTCUC 


J 


2.062 


1.611 


356 3 


-9 804 




2 


1.945 


1.537 


0.6485 


-18.1 




-J 


9 648 


9 930 


166 3 


1 859 


SFlora 


1 


2.089 


1.695 


345.6 


-3.195 


SFlora 


2 


1.865 


1.416 


349 


-7.847 


SFlors, 


3 


2.363 


1.986 


186.3 


5.415 


QMpH s 


1 


2.448 


2.017 


207.6 


6.152 


9Metis 


2 


2.619 


2.42 


199.1 


2.71 


QMetis 


3 


2.635 


2.463 


305.8 


-3.366 


9Metis 


4 


2.468 


2.067 


298.4 


-6.367 




1 


3.245 


3.09 


141.5 


-1.422 


1 OT-rvQipa 


2 


3.07 


2.675 


132 


-3.253 


1 OT-Tvcipn 


3 


2.795 


2.455 


235.1 


-4.039 


1 Kjny gicd 


4 


z,. / OJ 


9 4SS 


990 5 


-2.321 


1 1 Partlipnonp 




2.531 


2.303 


98.14 


-3.939 


1 1 ParfViPTior^p 


2 


2.658 


2.316 


88.96 


-1.312 


1 1 Par+ViPtinnp 


3 


2.633 


2.222 


200.7 


4.421 


1 1 PnT+ViPtinnp 

L ±1 ill tlldlU L/C 


4 


2.489 


2.255 


192.2 


5.076 


1 9Virtnria 




2.299 


1.906 


232.1 


-4.832 


1 9Virtnria 


2 


1.895 


1.585 


233.6 


5.213 


1 9Virtr»ria 


3 


2.386 


2.181 


78.76 


0.36 


1 9Virtnnfl 


4 


2.412 


9 066 


80 05 


-0 00617 


1 '^PcTPfia 


1 


2.686 


2.452 


269.3 


-7.478 


1 ^Ptjpna 


2 


2.781 


2.525 


260.9 


-15.85 


1 ^P<ypna 


3 


2.725 


2.3 


10.81 


-16.19 


13Egeria 


4 


2.598 


2.419 


1.3 


-6.083 


ISEunomia 


1 


2.928 


2.737 


278.5 


-7.263 


ISEunomia 


2 


2.625 


2.3 


269.6 


-0.2105 


15Eunomia 


3 


2.146 


1.877 


63.69 


13.24 


16Psyche 


1 


2.616 


2.197 


309.1 


-0.06303 


16Psyche 


2 


2.59 


2.309 


79.71 


-3.43 


16Psyche 


3 


2.785 


2.534 


74.13 


-2.755 


16Psyche 


4 


3.167 


2.794 


168.9 


0.01748 



Wyatt, M. C, Clarke, C. J., & Booth, M. 2011, Celestial Mechanics and 

Dynamical Astronomy, 111, 1 
Zacchei, A., Maino, D., Baccigalupi, C, et al. 2011, A&A, 536, A5 

' APC, AstroParticule et Cosmologie, Universite Paris Diderot, 
CNRS/IN2P3, CEA/lrfu, Observatoire de Paris, Sorbonne Paris 
Cite, 10, rue Alice Domon et Leonie Duquet, 75205 Paris Cedex 



18 



13, France 
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Planck Collaboration: Zodiacal Emission 
Table A.3. Parameters for observation seasons of Solar System objects detected by Planck. 



Object 


Season 


Solar Range 


Planck Range 


Ecliptic Longitude 


Ecliptic Latitude 






[AU] 


[AU] 


[deg] 


[deg] 


17Thptis 

1/1 llClla 


\ 


2.707 


2.379 


63.95 


-6.276 


1 / liiClla 


2 




2.62 


^'K 71 


-4.611 


1 / IIICUIS 


3 


2.643 


2.394 


155.5 


0.7942 


1/1 IlCllO 


4 




1 087 


145 6 


4 741 


1 SA^plnoTTiPTip 

J. OiVldlJlJllldlC 


\ 


1.894 


1.632 


31.15 


-1 1.75 


1 SlVTplnoTTiPTip 

J. OiVlCllJlJliidlt- 


2 


2.66 


2.279 


184.3 


2.536 


1 SlVfplnompiiP 

J. OiVlV^llJlJliit'llt' 


3 


2.792 


2.47 


169.8 


7.483 


1 QPnrtiiTia 


1 


2.067 


1.709 


77.41 


-0.4981 


1 9Pnrtiitifl 


2 


2.263 


1.938 


nil 


-1.697 


1 0T^nrtiina 

1 jfi- yji. uuiici 


3 


2.757 


2.354 


203.5 


-0.8577 


1 QPnfhina 

1 ^M. Kfl. UUllCl 


4 


2.828 


2.612 


192 


03586 




1 


2.301 


1.983 


357.8 


-0 01093 


i-V/-LV±Cli3i3 tux CI 


2 


2.154 


1.743 


177.3 


-0.7655 




3 


2.384 


2.015 


171.7 


-0.1937 


'70 A rnnVi i tntp 
^jy/T-iiiL/iiiii lie 


\ 


2.725 


2.532 


283.5 


-6.528 


9Q Amnhitritp 


2 


2.638 


2.3 


276.7 


-5.861 


90 Amnhitritp 


3 


2.406 


2.121 


44.22 


2.694 


41 T^nnlinp 

1 L/dJJIIIlC 


J 


3.226 


2.903 


345.7 


-1.722 


41 Dflnhnp 

r J- j^c^L/J-iiiw 


2 


3.516 


3.223 


56.65 


-11.57 


41 F)ar»linp 


3 


3.461 


3.335 


43.09 


-14.57 


41 Daphne 


4 


2.986 


2.639 


122.8 


-17.39 


4^PnfTpni n 

1 ^ U. tik-'llitl 


1 


2.551 


2.362 


300.8 


5.514 


4S Pi I OPT! 1 a 

1j LlgCl lid 


2 


2.679 


2.279 


296.1 


1.392 


4Sl-' 11 fTj^Ti 1 n 
J 1j LlgCllld 


■X 






41.13 


-■S 7 


4S Pi I OPT! 1 n 

1j LlgCllld 


4 


2.945 


2.607 


30.85 


-7.412 


S9Piiror>a 




2.847 


2.62 


88.9 


-7.082 


^^Piirnnfl 


2 


2.771 


2.462 


82.76 


-3.525 




3 


2.906 


2.525 


191.7 


5.929 




4 


3.065 


2.841 


183.3 


7.741 


OO 1 lllOUC 




2.391 


1.964 


329.2 


6.137 


OO 1 lll&UC 


2 


2.901 


2.667 


89.38 


2.673 


RSThisbe 


3 


3.108 


2.838 


78.56 


-0.01454 


RSThisbe 


4 


3.216 


2.981 


159.3 


-3.99 


88Thisbe 


5 


3.213 


2.828 


160.4 


-4.316 


128Neiiiesis 


1 


2.493 


2.23 


89.49 


-1.293 


1 ^R^PTYIPQIQ 
IZ.Ol'l ClllCOlo 


2 


2.683 


2.356 


82.97 


3.426 


lZ,Oi> ClllCdlO 






9 671 


1 84 'i 

lot. J 


7 flS4 


1 ^RN^pmpGiQ 

1^0i-> &111&D10 


4 


3.091 


2.801 


172.6 


6.189 


JZ.'-rlJ dlllUC'l gd 




1.993 


1.675 


97.96 


13.12 


JZiT-lJ dlllUC'l gd 


2 


2.572 


2.204 


92.5 


7.127 


J zJ-rQ dlllUCl gd 


■I 


'^08 




181 

101 


59^ 


■J^T^lJdlllUdgd 


4 


3.566 


3 227 


165.4 


-6.94 


J 1 il^dViU-d 




3.422 


3.071 


998 5 


17.37 


^ 1 ll-'dVllJ.d 


2 


3.626 


3.496 


217.5 


13.4 


O 1 il-'dViU.d 




D, I o\y 


^ SOS 


980 ^ 


6 678 


'S 1 1 DaviHa 

O 1 IL'dVlUd 


4 


3 717 


3.349 


268 


1.344 


704Ttifprannni a 

/ »jT"lIIld dlllllld 




3.498 


3.201 


230.6 


-17.46 


704Interamma 


2 


3.381 


3.248 


219.7 


-12.38 


704Interamnia 


3 


3.053 


2.907 


297.9 


-0.6138 


704Interaninia 


4 


2.806 


2.43 


289.8 


10.38 


Christensen 


1 


3.235 


2.864 


290.7 


22.74 


Christensen 


2 


4.137 


4.038 


300.2 


-10.42 


Christensen 


3 


5.066 


4.911 


261.8 


-24.26 


Christensen 


4 


6.385 


6.318 


273.7 


-34.56 


Christensen 


5 


7.408 


7.222 


245.3 


-40.48 



Westville Campus, Private Bag X54001, Durban 4000, South Africa 

Atacama Large Millimeter/submillimeter Array, ALMA Santiago 
Central Offices, Alonso de Cordova 3107, Vitacura, Casilla 763 



0355, Santiago, Chile 

' CITA, University of Toronto, 60 St. George St., Toronto, ON M5S 
3H8, Canada 

CNRS, IRAP, 9 Av. colonel Roche, BP 44346, F-31028 Toulouse 
cedex 4, France 

California Institute of Technology, Pasadena, California, U.S.A. 

Centre for Theoretical Cosmology, DAMTP, University 
Cambridge, Wilberforce Road, Cambridge CB3 OWA U.K. 

Centro de Estudios de Fisica del Cosmos de Aragon (CEFCA), 
Plaza San Juan, 1, planta 2, E-44001, Teruel, Spain 

14 r"^,,.,, — 1 r^^^™^i — — t- — t — , r> — — . m^^: 1 
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Table A.4. Flux density in Jy (IRAS -convention) of Solar System objects detected at a significance greater than 3o-, using PSF 
fitting. Error bars come from the statistical errors of the PSF fit and a calibration error of 1%, 5%, and 5% at 353, 545, and 857 GHz, 
respectively. 



Object Frequency Season 1 Season 2 Season 3 Season 4 Season 5 



1 Ceres 


857 GHz 




1 Ceres 


545 GHz 




1 Ceres 


353 GHz 


2.6 ± 0.48 


2Pallas 


857 GHz 


8.9 ± 0.076 


2Pallas 


545 GHz 


4.31 ± 0.11 


2Pallas 


353 GHz 


1.7 ± 0.19 


3Juno 


857 GHz 


2.96 ± 0.074 


BJuno 


545 GHz 


1.13 ± 0.36 


4Vesta 


857 GHz 


7.39 ± 0.068 


4Vesta 


545 GHz 


2.93 ± 0.1 


4Vesta 


353 GHz 


1.3 ± 0.24 


6Hebe 


857 GHz 


1.9 ±0.1 


SFlora 


857 GHz 


1.19 ± 0.11 


9Metis 


857 GHz 


1.19 ±0.19 


lOHygiea 


857 GHz 


2.5 ±0.13 


lOHygiea 


545 GHz 


1.18 ±0.17 


12Victoria 


857 GHz 




13Egeria 


857 GHz 




ISEunomia 


545 GHz 




16Psyche 


857 GHz 


1.18 ± 0.38 


ISMelpomene 


857 GHz 


1.05 ± 0.27 


19Fortuna 


857 GHz 


2.63 ±0.11 


20Massalia 


857 GHz 




29Amphitrite 


857 GHz 




45Eugenia 


857 GHz 


1.54 ±0.28 


52Europa 


857 GHz 


2.81 ± 0.12 


52Europa 


545 GHz 


1.08 ± 0.26 


88Thisbe 


857 GHz 


1.59 ± 0.33 


128Nemesis 


857 GHz 


1.4 ±0.27 


324Bamberga 


857 GHz 


2.83 ± 0.1 


SllDavida 


857 GHz 




704Interamnia 


857 GHz 


1.61 ± 0.29 


Christensen 


857 GHz 


2.6 ± 0.19 



17.1 


± 0.096 


17.7 ± 0.07 


16.1 ± 0.086 


6.72 


± 0.11 


7.51 ± 0.084 


6.68 ± 0.11 


2.52 


±0.15 


2.99 ± 0.094 


2.5 ± 0.14 


5.08 


± 0.072 


3.05 ± 0.094 


4.2 ± 0.13 


2.27 


±0.15 






1.79 


±0.23 


0.99 ± 0.13 




9.32 


+ 0.066 


10.1 ± 0.058 


8.29 ± 0.055 


3.55 


± 0.087 


3.96 ± 0.086 


3.18 ± 0.083 


1.58 


±0.32 


1.8 ±0.15 


1.2 ± 0.19 


1.52 


± 0.085 






1.78 


+ 0.089 












1.06 ± 0.24 


4.03 


± 0.11 


4.5 ± 0.1 


4.56 ± 0.093 


1.61 


± 0.24 


1.81 ± 0.25 


2.11 ± 0.15 


1.11 


±0.23 












1.17 ±0.16 






1.7 ± 0.33 




1.47 


± 0.21 


1.12 ±0.2 




0.93 


± 0.22 










0.913 ± 0.17 




1.21 


±0.19 






3.38 


±0.26 


2.05 ± 0.12 


1.57 ±0.11 


1.79 


± 0.13 










1.07 ± 0.32 




1.42 


±0.17 


2.19 ±0.16 


1.88 ±0.15 
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Fig.A.2. Top Panel: the ratio of / ■ ^fs of the first and subsequent measurements for those Solar System objects which were 
detected in more than one survey. Measurements at 857 GHz are shown in blue, measurements at 545 GHz are shown in green, 
and measurements at 353 GHz are shown in red. The grey line shows (d/do)^, which one would expect for Rayleigh-Jeans ob- 
jects. Bottom Panel: the spectral index for those asteroids detected in multiple frequency bands for a given survey, defined as 
log (7545/7857) / log (545/857). Juno and Vesta are shown in brown and gold, to indicate that they have different spectral classifi- 
cations than Ceres, Pallas, Hygiea and Europa, shown in shades of blue. The white symbols show the corresponding values for 
353-545 GHz, when they exist. For both the top and the bottom panels, the unique symbols used for each object are: 1 Ceres - 
circles; 2 Pallas - squares; 3 Juno - diamonds; 4 Vesta - hexagon; 6 Hebe - -n; 8 Flora - X; 9 Metis - upward-pointing triangle; 
10 Hygiea - right-pointing triangle; 19 Fortuna - vertical line; 45 Eugenia - downward-pointing triangle; 52 Europa - star; 324 
Bamberga - left-pointing triangle; 704 Interamnia - horizontal line 
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